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Abstract: The MEMS vector hydrophone is a novel acoustic sensor with a "four-beam- 
cilia" structure. Based on the MEMS vector hydrophone with this structure, the paper 
studies the method of estimated direction of arrival (DOA). According to various research 
papers, many algorithms can be applied to vector hydrophones. The beam-forming 
approach and bar graph approach are described in detail. Laboratory tests by means of the 
a standing- wave tube are performed to validate the theoretical results. Both the theoretical 
analysis and the results of tests prove that the proposed MEMS vector hydrophone 
possesses the desired directional function. 

Keywords: MEMS; piezoresistive effect; vector hydrophone; DOA estimation; 
beam-forming 
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1. Introduction 

Vector hydrophones can obtain both the acoustic pressure and acoustic particle velocity 
simultaneously, and therefore can gather more information about a point in the sound field than an 
ordinary acoustic pressure hydrophones. Thus, they have a unique advantage and wide applicability [1-4]. 

As early as in the 1940s, the United States developed a sound pressure gradient vector hydrophone, 
and in 1970s, the vector hydrophone was successfully applied in the DIFAR sonobuoy system. In the 
development of vector hydrophones, the United States and Russia have taken the lead in the World, 
and vector hydrophones with stable performance have not only gone into the engineering phase [5-6], 
but explorations in the field of vector hydrophone calibration have also been carried out [7]. In the 
field of the vector hydrophone signal processing, Nehorai established a position measurement model 
with a combination of sensor arrays in 1994, and gave two methods to measure the orientation error, 
the methods are the Cramer-Rao lower limit and the use of a single combination of sensors, 
respectively [8]. Because traditional hydrophones are piezoelectric ceramics, their structure and 
principle are also relatively simple, so the United States and Russia also typically use the resonant 
vector hydrophone based on the piezoelectric accelerometer. 

The MEMS vector hydrophone is designed based on the piezoresistive effect of silicon, The beam 
microstructure is manufactured by means of a silicon-on-insulator (SOI) wafer with MEMS 
technology. Compared with traditional piezoelectric hydrophones, the MEMS vector hydrophone has 
the characteristic of a being a single sensor with directional functions. In theory, two sensors can locate 
the target. 

2. MEMS Vector Hydrophone 

In this paper, we will discuss a novel hydrophone which is developed from a typical accelerometer 
structure. The SEM image (top view) of the microstructure is shown in Figure 1. 

Figure 1. SEM images of the microstructure. 




As we can see in this Figure, it is a typical mass-beam accelerometer structure, which is based on 
the piezoresistive effect. In order to improve the sensitivity, a cylinder (diameter: 200 |um, 
length: 7,000 |um) was glued in the center of the mass [9]. Figure 2 shows the bionic microstructure of 
MEMS vector hydrophone, which consists of two parts: four high-precision cantilever beams and a 
rigid cylinder which is fixed at the central of the square mass. Figure 3 shows the frequency response 
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of this geometry from 0 Hz to 1 KHz, BK 8305 as a standard accelerometer (sensitivity is 60 mV/g), 
the sensitivity of this structure in the X direction is 0.7552 mv/g. 

According to the auditory principle of a fish's lateral line, we use a rigid cylinder as a stereocilia 
which can improve the sensitivity. When the sound signal is sensed by the cylinder, the piezoresistors 
located at the beam transform the sound signal into strain and finally into a differential output voltage 
signal via the Wheatstone bridge circuit. So the sound signal can be detected. 

Figure 2. Schematic view of the bionic structure of the MEMS hydrophone. 




Figure 3. Frequency response of the structure. 
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3. The Principle of Direction 

The sound field is the vector field, and a plane wave is a longitudinal wave, so the vector 
hydrophone measures the direction of sound vibration (vibration velocity v) and sound intensity flow P 
v (where P is acoustic pressure) is the direction of sound propagation, which is the goal orientation, so 
the measurement Pv can estimate the DOA. 

From the acoustic theory, we can see the plane wave acoustic pressure can be expressed as [10-12]: 

p(r,t) = p 0 e J{kr ~ wt) = p 0 exp[j (he cos a cos 0 + ky cos a sin 6 + kz sin a -wt)] (1) 
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Where k is a wave vector, indicating the direction of sound propagation, a is the angle between k 
and horizontal plane, a G [-nil, n/2], 9 is the angle between the projection in the horizontal plane of k 
and x-axis, 9 G [0,27i], As shown in Figure 4: 

Figure 4. Wave vector's projection in the coordinate system. 

z 




In homogeneous medium, the sound field equation of motion is: 



£ + -V,--0 (2) 
dt p 0 



Substituting this into Expression l,we get: 

v= — (cos^cosa^+sin^cosa/7+sinar) (3) 

Po c 

Where £ //, r is the unit vector of x, y, z axis, /?o is the medium density, c is the velocity of sound in 
the medium. Equation (3) shows that the ratio between acoustic pressure and the three components of 
particle velocity is just a constant, their wave is the same, so pressure and velocity of the plane wave 
are totally relevant. 

From (3) we can get the three components of the particle velocity: 

v x = — ^-cos#cosa 

p 

v = sin#cosa (4) 

P ■ 

v z =- £ — sin a 

and: 

0 = arctan(^) (5) 
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a = arctan(— = z ) 




(6) 



So, as long as we measured the two components v x ,v y of particle velocity in the horizontal plane, 
we can get the azimuth 6 of the sound source in the horizontal plane from Equation (5), three velocity 
components can get the pitch angle a from Equation (6). These are the basic principles that the vector 
hydrophone uses to determine a sound source's position. 

4. Two Approaches to DOA Estimation 

According to the different types of noise background and signals, there are several approaches to 
estimate azimuth angle based on a single vector hydrophone [13], such as the average acoustic power 
approach, line-spectra estimation, the bar graph approach and the beam-forming approach. However, 
the first two methods can't do well in the condition of either line-spectra coherent interference lying in 
broad band signal or broadband coherent interference lying in line-spectra signal. In this paper, we 
apply the beam-forming approach and bar graph approach to the DOA estimation based on a single 
MEMS vector hydrophone. 

4.1. Beamforming Approach 

The term beam-forming refers to the various output signals whereby a vector hydrophone can 
determine space directivity after processing (weighted, delay, summation etc.), It can focus the receive 
direction of hydrophone on a direction, which is equivalent to a beam, rotating the beam to find the 
maximum peak, that is the target azimuth. The algorithm can be seen as a space filter that can separate 
the signal from interference, according to the azimuth difference between signal and interference. 

There are many weighted modes, this paper only analyzes two of them. If only in the plane, 
according to (4) we can obtain the following equation: 




(7) 



P 



sin# 



The first weighted method: v x (t) cos y/ + v y (t) sin y/ 

Substituting (7) into it and assuming — !— =1, we have the following equation: 

v x (t) cos y/ + v (t) sin y/ = p x cos(# - y/) 



(8) 



The second weighted method, adds the acoustic pressure p(i), so we get: 

P + v x (0 cos ¥ + v y (0 s i n ¥ = PQ + cos(# - ¥)) 



(9) 



The directivity patterns of two weighted functions (8) and (9) are as follows, respectively: 
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Figure 5. Directivity patterns of function (8) (6 = 0°and 6 = 30°). 




The Figure 5 is "8"-shaped with bilateral directivity, which can induce port and starboard ambiguity, 
however, the Figure 6 is heart-shaped with unilateral directivity. Therefore, we use the second 
weighted form to reduce the side-lobe. If we don't add p(i), the outcome will have two peaks. 

The principle chart of the beam-forming approach is shown in Figure 7, weighting the three way 
signal p{t), v x (t) and v y (t) output of hydrophone by 1, cos#, sin# and summing them. The result is: 

y(f) = p(t) + v x (t) cos 0 + v y (t) sin 0 (1 0) 



The average power of y (t) is: 



M{6) = Ey(t) 



(11) 



where, E[ ] — ensemble average; w(ff) — space spectrum of the output, is a function of 9 , reflecting the 
space distribution of energy. Search the max of w{6) in #^[0,27i], w{6) can acquire the max, only 
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when 6 = #o, so 9o is the azimuth angle of the sound signal. Virtually, this is a process of beam- 
forming [14]. 

Figure7. Principle chart of the beam- forming approach, 
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4.2. Bar Graph Approach 

From Equation (5) we know that, by simply calculating the three velocity components measured by 
a vector hydrophone, we can identify the azimuth angle 6 and pitch angle a of a sound source. By 
determining the positive and negative of v x and v y , a single vector hydrophone can estimate the 
direction of the sound source. Although this method is simple, the error would be great, this is mainly 
because the three velocity components v x ,v y and v z , measured by the vector hydrophone not only 
contain a valid signal, but also noise. Noise has a great impact on the estimated results, and this impact 
is more serious with lower SNR (signal-to-noise ratio). 

We assume the signals received from vector hydrophone are as follows: 

p(t) = p s (t) + p n (t) (12) 
v(t) = v s (t) + v n (t) (13) 

In the equation, p s (t)- acoustic pressure generate by target sound source; 

v s (0- particle velocity generate by target sound source; 

p n (t)- acoustic pressure of isotropic noise field; 

v n (t)- particle velocity of isotropic noise field. 
Target signal and interference noise is independent of each other, and they are ergodic. Calculating 
average sound intensity by above formula, we get: 

W) = Pio ■ m = (p s (o + p n (0) • (v, (o + v n (0) * Ps (ov, (o (H) 

For p n (t), v n (t) and p s (t) are independent of each other, the value of the average cross-multiplied is 
small, which can be neglected. 



0 = arctan 



hit) 



(15) 



Bar graph estimation is a statistical approach based on cross-spectrum estimation. First we calculate 
the angle corresponding to each frequency, and then count the probability density of each estimated 
value, to get the estimated curve of a certain moment, the curve maximum position corresponds to the 
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estimated value of which is the target position. Figure 8 shows its principle chart [13], acoustic 
pressure and particle velocity are calculated by the conjugate and cross-spectrum. If P(f)> V x (f) and 
Vyif) correspond to the FFT of the acoustic pressure and particle velocity signal p(t)^ v x (t) and v y (t), 
we can get the expression of calculated azimuth: 



6 = arctan 



I y(t) 



= arctan 



R[P*(f)v y (f)] 
R[P*(fW x (f)] 



(16) 



where, R[ ] — get real component; "*" — conjugate. 

Bar graph statistics mainly estimates the azimuth of every frequency, if 1° is the statistical unit, we 
get: 

k = [0(f)xlSO/7r] (17) 
<p(k) = <P(k) + \ (18) 

In Formula (17), [ ] — get integer. 9(f) is the azimuth value of certain frequence, transforming 9(f) 
to angle k. cp is an array, which stores the frequence of every angle in [-180°, 180°], initialization of the 
array is zero, when we get a value of k, the frequence cp(k) of corresponding to k add one in array. The 

statistical result is the azimuth estimated value. 

Figure 8. Principle chart of the bar graph approach. 
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5. Test and Results 

The DOA estimation experiment has been completed by means of the a test vector hydrophone 
instrument in the laboratory, Figure 9 shows the photo of MEMS vector hydrophone being tested. Its 
receiving sensitivity is up to -160 dB(0dB = 1 V/|upa). Figure 10 shows the principle chart of the test 
instrument, the instrument can emit standing wave in a circular tube. The sound source is located in the 
tube's bottom, so a uniform and stable sound field is built in the tube. We can change the angle of 
wave incidence by adjusting hydrophone's position. The computer collects the sound information 
transformed by hydrophone, and processes the data, and as a result, the direction of the sound signal 
can be estimated. 



Sensors 2009, 9 



6831 



Figure 9. The photo of the MEMS vector hydrophone. 




Figure 10. Principle chart of the calibration instrument. 
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5.1. Directed to the Single-frequency Signal 

Using a signal generator to generate a single-frequency sine wave Si(t), DOA is 45°. The Data 
Acquisition Card collects a three-way signal at the same time P, V x , V y , data sampling rate is set 
to 10 KHz. The two methods are applied to DOA estimates. The results in Figure 11 show Method 1 
for the beam-forming approach, Method 2 for the bar graph approach. 
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Figure 11. The results of the two methods (DOA = 45°). 
Method 1: angle = 46.62° Method 2: angle = 46° 
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5.2. Directed to the Noise Signal 

As sound source a broadband signal (white noise) S2OO is used, DOA is 0°. The results of the two 
methods are shown in Figure 12: 



Figure 12. The results of the two methods (DOA=0°) 
Method 1: angle = 0.18° Method 2: angle = -1° 
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5.3. Directed to the Line-Spectrum Signal in the Context of Noise 

As sound source we enter a single-frequency signal with attached white noise Si(t) + S2C:), DOA 
is 90°. The results of the two methods are shown in Figure 13: 
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Figure 13. The results of the two methods (DO A = 90°). 
Method 1: angle = 100.98° Method 2: angle = 89° 
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6. Conclusions 

In this paper, by testing and analyzing the directional functions of a MEMS vector hydrophone in an 
ideal room environment, we obtained better results. From the test results, the MEMS hydrophone can 
not only direct but also have high orientation accuracy. For a continuous wave, two methods are 
feasible, but because of the requirements of frequency resolution, the bar graph approach is not 
suitable for pulse waves. There are two main reasons of error: first, because of the inaccuracy of the 
measured angle and the other is the rounding error of the bar graph. It remains to be determined 
whether the MEMS vector hydrophone can achieve good goal orientation or not in the actual marine 
environment. 
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